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We present a theoretical analysis of the effect of inelastic electron scattering on current and its 
fluctuations in a mesoscopic quantum dot (QD) connected to two leads, based on a recently developed 
nonperturbative technique involving the approximate mapping of the many-body electron-phonon 
coupling problem onto a multichannel single-electron scattering problem. In this, we apply the 
Biittiker scattering theory of shot noise for a two-terminal mesoscopic device to the multichannel 
case with differing weight factors and examine zero-frequency shot noise for two special cases; (i) 
a single-molecule QD and (ii) coupled semiconductor QDs. The nonequilibrium Green's function 
method facilitates calculation of single-electron transmission and reflection amplitudes for inelastic 
processes under nonequilibrium conditions in the mapping model. For the single-molecule QD we 
find that, in the presence of the electron-phonon interaction, both differential conductance and 
differential shot noise display additional peaks as bias-voltage increases due to phonon-assisted 
processes. In the case of coupled QDs, our nonperturbative calculations account for the electron- 
phonon interaction on an equal footing with couplings to the leads, as well as the coupling between 
the two dots. Our results exhibit oscillations in both the current and shot noise as functions of 
the energy difference between the two QDs, resulting from the spontaneous emission of phonons in 
the nonlinear transport process. In the "zero-phonon" resonant tunneling regime, the shot noise 
exhibits a double peak, while in the "one-phonon" region, only a single peak appears. 

PACS numbers: 73.63.Kv, 71.38.-k, 73.50.Td 



I. INTRODUCTION 

Recent progress in nanotechnology has made possible 
the fabrication of single-electron tunneling devices us- 
ing organic molecules, the size of which is sufficiently 
small so that the discrete nature of its energy levels is 
important. In contrast to semiconductor quantum dots 
(QD), the molecular materials possess much weaker elas- 
tic parameters, such that it is very easy to excite their 
internal vibrational degrees of freedom (phonon modes) 
when electrons are incident upon the molecule through 
a tunnel junction from an external environment. This 
can significantly influence electron tunneling at low tem- 
perature, and the phenomenon has now provoked a large 
amount of experimental^'^''^'^'^ and theoretical^'^'^'^'^°'^^ 
work on the problem of tunneling through a single level 
with strong coupling to local phonon modes. 

On the other hand, phonon-asistcd inelastic tunneling 
in semiconductor QDs has continued to attract extensive 
studies. In particular, a recent experiment has measured 
the nonlinear tunneling current through a tunable two- 
level quantum system, i.e., coupled semiconductor QDs 
(CQDs), at low temperature with the observation of ro- 
bust spontaneous emission of phonons leading to an os- 
cillatory structure in the current spectrum. Following 
this, nonperturbative theoretical analyses ascribed this 
oscillation to an interference effect of the electron-phonon 
interaction (EPI) in CQDs.^^^w 

Thus far, most theoretical work concerning this inelas- 
tic tunneling problem has focused on calculating cur- 
rent through a two-terminal setup in dc and ac con- 



ditions. Only a few attempts have been made to 
study the current-current correlations (shot noise) in QD 
systems, which has become a central subfield of 
mesoscopic physics. ^^'^^'^^ In the present paper, we in- 
vestigate low-frequency shot noise behavior when the 
QD is subject to strong coupling with a local phonon 
mode. This kind of strongly nonequilibrium problem has 
been found to be difficult to solve because previous stud- 
ies demonstrated that when inelastic processes such as 
phonon emission and absorption are considered in the 
tunneling event through a QD, electron scattering be- 
comes a strongly correlated many-body problem involv- 
ing electrons and the various excited phonon states of the 
QJ3 6,7,8,9,10,11,13,14 ^Yiis case, the usual perturbation 
theory is invalid for dealing with this scattering prob- 
lem, albeit that it has been extensively exploited with 
enormous success to study conventional transport in bulk 
materials. 

Recently, a nonperturbative scheme was proposed by 
Bonca and Trugman^^ to study the small polaron of the 
coupled electron-phonon system described by the Hol- 
stein model. Later, this method was applied to the anal- 
ysis of inelastic electron scattering in mesoscopic quan- 
tum transport in semiconductor QDs and molecular wires 
within the framework of the Landauer-Biittiker scatter- 
ing theory.^^'^^'^^ Very recently, we advantageously em- 
ployed this mapping scheme to evaluate time-dependent 
phonon-assisted tunneling at low temperature through 
a single-molecular QD with the help of the psuedo- 
operator technique and the nonequilibrium Green's func- 
tion (NGF) method.^'' The fundamental idea is to re- 
model the strongly correlated Hamiltonian in terms of 
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the combined electron-phonon Fock space, mapping the 
many-body problem onto a one-body scattering problem 
(noninteracting system), which is highly desirable to fa- 
cilitate an accurate calculation of the current-fluctuation 
spectrum. 

However, this mapping technique introduces an addi- 
tional difficulty: it also maps the ordinary single-mode 
leads in the two terminals onto pseudo-multi-channel 
leads associated with differing weight factors P„ cor- 
responding to the statistical probability of the excited 
phonon number state n of the mesoscopic device. Con- 
sequently, the original Landauer current formula needs 
to be modified to properly account for this fact and the 
Pauli exclusion principle as it applies to multi-channel 
leads. ^^'^"^'^^ For example, if we consider a tunneling 
event in which an electron with energy e enters from the 
nth channel of the left lead, scatters with a phonon in the 
QD, and finally exits into the mth channel of the right 
lead with energy e', the contribution of this process to 
current is simply proportional to P„/£(e)[l — (e')] 
is the Fermi distribution function of the Ith pseudochan- 
nel in lead rj, as described in Section II below), multiply- 
ing the square modulus of the transmission amplitude, 
liLri.flmP- Correspondingly, it is clear that an analogous 
modification must be made in the Biittiker scattering for- 
mula for current-current corrclations^^'^^ in the presence 
of multichannel leads with differing weight factors. This 
issue is a central focus of this paper. In Section II. we will 
briefly present the appropriate generalization of scatter- 
ing theory and provide the general formula for determin- 
ing the zero- frequency shot noise spectrum of mesoscopic 
systems with EPI. 

In two applications of this generalization we will ex- 
amine low frequency electric current fluctuations as func- 
tions of voltage for the two cases mentioned above: 1) a 
single-molecular QD (single site) and 2) semiconductor 
CQDs (two sites) in the central region sandwiched be- 
tween two normal leads. Every QD is taken to have only 
one single-particle energy level coupled with a disper- 
sionless phonon mode of the Holstein-type model. Here 
we concentrate our attention on the effect of the inelas- 
tic scattering process in tunneling and ignore all other 
complexity of real molecular or semiconductor devices, 
except for EPI. For the single site case, our calculations 
predict (i) a step-type characteristic of the bias voltage- 
dependent shot noise analogous to the current-voltage 
curve at low temperature, in which the first step cor- 
responds to resonance with the single level of the QD 
and the others result from phonon-emission-assisted res- 
onances, and (ii) an enhancement of the Fano factor due 
to phonon effects in the high bias-voltage region. 

For the two-site case, in order to observe the effect of 
spontaneous phonon emission, we consider an extraordi- 
narily high bias- voltage between the two leads and calcu- 
late the current and shot noise as functions of the energy 
difference between the two QDs, which can be tuned ex- 
perimentally using gate voltages. Very recently, the shot 
noise spectrum in CQDs was explored by means of master 



equations with correlation functions in Laplace transform 
space. Under the physical presumption of weak tunnel- 
ing between the QDs and leads, this method employed 
the noninteracting blip approximation to deal with the 
EPI.^'^ Our objection here is to improve upon this by 
developing an approach which allows us to account for 
the EPI on an equal footing with the couplings to the 
leads, as well as treating the coupling between the two 
dots in a nonperturbative way. The mapping technique 
conforms to this purpose. Our calculations exhibit an 
oscillatory shoulder in the current and shot noise spectra 
due to spontaneous phonon emission in nonlinear trans- 
port. In the resonant tunneling regime, the shot noise 
has a double peak, while only a single peak shows up in 
the phonon-assisted resonant tunneling region. 

The rest of this paper is arranged as follows. In Sec. 
II, we describe the class of systems that we study and 
present the appropriate generalization of scattering the- 
ory for calculating current fluctuations of a two-terminal 
device in the presence of EPI. In general, we find that 
the current fluctuation spectrum cannot be expressed in 
terms of transmission and reflection probabilities even 
in the zero frequency limit, but it is sensitive to the 
transmission and reflection amplitudes involving all pseu- 
dochannels. Our two applications occupy the following 
two sections: Sec. HI for a single molecular QD and 
Sec. IV for the CQDs. In these two special cases, we 
employ the NGF technique to evaluate transmission and 
reflection amplitudes in the wide band limit. Numerical 
results and some discussions are also presented in these 
sections. Finally, our conclusions are summarized in Sec. 
V. 



II. METHOD AND FORMULATION 

A. Model and mapping Hamiltonian 

Here, we consider the model of the two-terminal meso- 
scopic device as constituted of a central region (the de- 
vice) connected to two noninteracting reservoirs via tun- 
neling and coupled to a phonon bath. We assume that the 
two leads remain in local equilibrium with the Fermi dis- 
tribution, fr,{uj) = [l-fe^'^"'^'')/'^^'^]"^ [^ri is the chemical 
potential of the lead r] {= L or R), T is the temperature], 
and both of them have only one electronic channel. We 
also assume that the EPI takes place solely in the device 
and not in the ideal leads, specifically, there is no ex- 
change between electrons and phonons in the tunneling 
processes. This is reasonable because high-order tunnel- 
ing processes accompanied by phonon emission and ab- 
sorption are much weaker than direct tunneling events. 
For simplicity, we adopt a single-mode Einstein-phonon 
with frequency uiph- Moreover, we neglect the spin degree 
of freedom and any effects of electron-electron Coulomb 
interactions. 

Therefore, the Hamiltonian of this system involves 
three terms: H = iJioad + H^cn + Ht\, where H\ca.d de- 
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scribes the two isolated leads, Hcen models the interact- 
ing central region in addition to the free-phonon contri- 
bution, and Hti is the tunnel coupling between the leads 
and the device, respectively: 



^ccn 



q ., k 

Hd{{d]}-{d^]) + kujphb''b 



H.c). 



(la) 



(lb) 



(Ic) 



Here cj^j. {crik) are the creation (annihilation) operators 
for the noninteracting electrons with momentum k and 
energy e^^ in the lead 77, and {d]} {{dj}) forms the com- 
plete and orthonormal set of single-electron creation (an- 
nihilation) operators in the central region, respectively. 
The form of Hd depends on geometry of the device being 
investigated. The last term in Eq. (lb) denotes the inter- 
action of an electron on the central site with phonons: 
(6) creates (destroys) a phonon, and Xj is the on-site EPI 
constant, l^.j stands for the tunnel coupling between the 
device and the lead 77. 

Obviously, the model described by the above Hamil- 
tonian (1) involves a many-body problem with phonon 
emission and absorption when the electron tunnels 
through the central region. Following the suggestion of 
Bonca and Trugman, we expand the electron states in 
the central region and the leads in terms of the direct 
product states composed of single-electron states and n- 
phonon Fock states: 



|j,n)=rft^^|0), (n>0) 



and 



\r]k,n) 



|0), 



(2) 



(3) 



such that the electron state |j) in the central region and 
the state \r]k) with momentum k in lead 77, respectively, 
are accompanied by n phonons (|0) is the vacuum state). 
After performing this transformation, the many-body on- 
site EPI in Eq. (lb) can be mapped onto a one-body 
model:i9^2i,22,23,24 



E 

j,n>Q 



-AjVn + l(|j,n + l>(j>| + |j,7i)(j,n + 1|)]. (4) 



In this procedure, the two noninteracting single-mode 
leads of Eq. (la) are mapped to a multichannel model 



H\cad = ^ er,kn\rik,n){iik,n\ 



(5) 



with €rikn = ^Tjk + nfiLOph- Here, the difference from 
traditional multi-channel leads lies in the fact that the 



channel index n represents the phonon quanta excited in 
the device, thus generating a weight factor P„ = (1 — 
g-/3wph^g-n/3wph j^jjg statistical probability of the phonon 
number state \n) at finite temperature T {[3 = l/fceT)], 
which must be associated with this channel. In this con- 
nection, we assume that the phonon bath is perpetually 
in thermal equilibrium, even under nonequilibrium trans- 
port conditions, due to a high energy relaxation rate. 
Thus, we ignore nonequilibrium phonon effects, consid- 
ering phonon equilibration fast compared to the dwell 
time of an electron on the QD.^^ We also note that dif- 
ferent channels are admixed with each other due to the 
EPI term, Eq. (4), which takes place solely in the cen- 
tral region. Finally, the tunneling part (Ic) can also be 
rewritten in terms of this basis set: 



^ti= E V:i^{\i^k,n){]M 

ri.k,j,n 



H.c. 



(6) 



is the coupling between the nth pseudochannel in 
lead 7] and the device. Clearly, the many-body EPI prob- 
lem is transformed to a multichannel single-electron scat- 
tering problem with the help of the new representation, 
as illustrated in Fig. 1. 
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FIG. 1: Schematic description of the inelastic scattering prob- 
lem for (a) a single site and (b) two sites in series, both with 
on-site EPI. Each phonon state of the site, taken jointly with 
the Bloch state of the electron in the lead, can be visualized 
as a pseudochannel labeled by n, which is connected to the 
two leads with hopping parameters Ti^/ji. These channels are 
connected vertically by the EPI, Xj. 

We depict the physical process of tunneling after re- 
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modeling the electron-phonon coupling system as follows: 
The nth pseudochannel in either the left or the right 
lead denotes an electron propagating in one of the single- 
electron modes of that lead with the device being in its 
n-phonon state. An electron incident on pseudochannel 
n in the left or right lead can transport clastically to 
the central region, where it then excites or absorbs m 
phonons and exits inelastically into the (n ± m)th chan- 
nel in the left or right lead; or experiences no exchange 
with phonon and exits elastically in the same channel 
in both leads. It is of central importance to note that, 
within this formulation, the many-body problem can be 
solved exactly and consequently there is no loss of phase 
coherence in the treatment of the electron-phonon ex- 
change. Indeed, such coherent interference effects in the 
emission processes of phonons play a crucial role in non- 
linear transport in CQDs. This is the reason that the 
present approach can provide a satisfactory description 
of the current spectrum consistent with actual experi- 
mental results. 



B. Current-current fluctuation spectra 

In the section, we derive a general formula describ- 
ing low-frequency current fluctuation spectra for a two- 
terminal mesoscopic device coupled to an Einstein- 
phonon bath (with frequency coph). In the combined 
electron-phonon representation, the many-body problem 
becomes a one-body multichannel scattering problem, as 
described in the last section. With this picture in mind, 
we employ the Biittiker scattering approach to achieve 
this objective. ^^'^^ 

First we introduce creation and annihilation operators 
of electrons for the incoming and outgoing states in all 
transverse channels. We employ creation, ajj„(w), and 
annihilation, a^„(a;), operators for electrons with total 
energy lo incident upon the central region from channel 
n of lead rj; and we employ the creation, ^„(cj), and an- 
nihilation, &^„(w), operators for electrons in the outgoing 
states. Clearly, they obey anticommutation relations. 

In second quantized notation, we can write explicit 
forms for the wave function \E',,(r,t) in lead 77 as: 



*r,(r, t) 



duje 



'iujt/h 



N 

E 



Xr,n(rj 



(7) 



Here, ^re the transverse parts of the electron wave 
hmctions [r = (rj^, z)], the wave vector in the nth chan- 
nel of lead 77 is = h~^\p,m* {uj — nhiOph)Y^'^ {'m* is 
the effective electron mass), and the corresponding elec- 
tron velocity is w^„(w) hk^pi/m* . N denotes the total 
phonon number, which is infinite in principle, but must 
truncated to a finite value for computational purposes. 
Correspondingly, we define the current operator as: 



- , , he 
Zim 



dr I 



dz 



(8) 



Substituting the wave function of Eq. (7) into this defini- 
tion, and assuming that the velocities w^„(a;) vary with 
energy quite slowly, the current operator can be simpli- 
fied as: 



I,f{t) = T E/ / ^'^d.Uj't 
n '' 

-^L(^)^rm(^')]- 



,i(w-w')t/?ir„t 



(9) 



In the scattering approach, the outgoing electron opera- 
tors, b, are related to the incoming operators, a, via the 
2N X 2N scattering matrix s: 



s = 



sl.l sl.r 
sr.l sji^R 



tR Tr 



(10) 



in which all blocks have N x N dimensions. An element 
Srin,rim of thc diagonal blocks describes a tunneling event 
in which an electron is incident from the 77m channel and 
reflects back to the nth channel in the same lead. For 
this reason we denote it as rLnm {rRnm), the reflection 
amplitude for the left (right) lead. Similarly, an element 
Snn,r)'m {v' 7^ ^7) of the off-diagonal blocks denotes elec- 
tron transmission through the central region from the 
?7'm channel to the 77n channel; and these are the trans- 
mission amplitudes t^nm and tRnm- It should be noted 
that the scattering matrix is unitary s^s = / (/ is the 
iV-dimensional unit matrix) and symmetric. Therefore, 
we can express the current in terms of the operators a 
and alone: 



E 

a.(3^mn 



xA-;(77;^,c.')a^„(c.'), (11) 



where 



(12) 

Proceeding to evaluation of the shot noise spectra, we 
recall that the noise power spectra between the leads 77 
and 77' can be expressed as the Fourier transform of the 
current-current correlation function: 

+ {ljit%{t))-2{i,{t)){i,,{t'))], (13) 

where (•••) represents the quantum statistical aver- 
age. Focusing attention on low-frequency noise in the 
same lead in the present paper, we obtain, after some 
algebra^^'^^ 

^.'^(0) = X E E / dujdu;'A::^{r^;cu,cu) 
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x4l(^;^',^'){(«L„(^)«w('-'))(a,J^)a;:,(c.')) 
(aU^>pS^))Ki{^>lM)}- (14) 



Considering the fact that the channel actually corre- 
sponds to the phonon number state in this mapping 
model, these quantum statistical averages, (a^a), must be 
calculated with caution: an additional weight factor P„ 
must accompany any starting channel n jointly with the 
Fermi distribution, /,^(w) = [1 + e('^+"''"'"'-'^'')/'==^^]~^ 
Therefore, we have 

X {P™/™(..)[1 - f^ico)] + P„/^(c.)[l - /™(^)]} . (15) 

In contrast to the case without EPI, one can not further 
simplify the shot noise formula of Eq. (15) to write it in 
terms of transmission probabilities only without any de- 
tailed information concerning the device. Here, the shot 
noise formula for inelastic tunneling requires information 
about the transmission and reflection amplitudes of all 
pseudochanncls. Actually, a similar result, that the noise 
can not be expressed in terms of the transmission prob- 
abilities only in the boson-field assisted tunneling, has 
been addressed in the studies of ac-driven current noise 
in quantum dots.^^ In the next two sections, we will ap- 
ply this theory to evaluate zero-frequency shot noise in a 
single-molecular QD and semiconductor CQDs based on 
Eqs. (12) and (15). 

It should be noted that since the mapping tech- 
nique involves an underlying single-particle picture for 
the dot EPI, it does ignore the possibility of strong 
bath-dot many-body effects which could invalidate the 
single-particle picture, introducing broadening of tun- 
neling associated with strong EPI and nonequilibrated 
phonon. ^'^^ Such many-body effects are understood to 
play an important role in tunneling in the linear regime, 
but become trivial in the case of high bias- voltages, upon 
which this paper is focused. 



III. A SINGLE-MOLECULAR QUANTUM DOT 

A. Model and shot noise formula 

In this section, we examine the EPI effects on the zero- 
frequency noise properties of a single-molecular QD. For 
this single-site case [shown in Fig. 1(a)], the Hamiltonian 
Hd of the central region becomes 



(16) 



and summation with respect to the site index j reduces to 
only one term. In terms of the combined clcctron-phonon 
representation, we rewrite it as: 



n>0 



nhiUph)\l,n){l,n\. 



(17) 



As suggested by Haule and Bonca in their original 
work,^^ the transmission and reflection amplitudes in 
this single-site case can be obtained by solving the 
Schrodinger equation with trial wave functions. In 
our recent study of low-temperature time-dependent 
phonon-assisted tunneling through a single-molecular 
QD, we adopted the NGF method to compute the time- 
depdcndent transmission probabilities within the frame- 
work of the mapping technique. In fact. Fisher and Lee 
established a general relation between the scattering ma- 
trix elements and the retarded GFs of the conductor. 
Hence, in the present paper we will employ the NGF 
technique to compute these transmission and reflection 
amplitudes in the wide band limit. 

Following our previous paper, ^'^ we define pseudo- 
Fermi operators: 



|l,n)(0|,c: 



7]kn 



\vk,n){0\. 



(18) 



and rewrite the mapping Hamiltonian in terms of these 
pseudoperators. The evaluation of the retarded GFs of 
the QD, G;„(i,i') = ie{t-t'){{di^{t), 4^{t')}), pro- 
ceeds using the equation-of-motion technique. In the 
Fourier space, we have 

[u -€d- mhhjph - i:',„(a;)] G'^„(a ;) = S^n 
-AiV^G'[„_i)„(w) - AiVm+ lG[,„+i)„(w),(19) 

where the retarded self-energy Y^I-^ {lu) due to coupling to 
the leads is defined as 



I ^rjm I 



^ ~ ^'i-jk ~ ITihuJpii + zO+ 



(20) 



In the wide band approximation, where the hopping ma- 
trix element V^^" ^ Vrij is independent of energy, one 
has = -i(r^ + Tr) with r„ = 27tJ2, \V„\^S{u; - 
^rjk — mhiOph) as the generalized linewidth function. For 
simplicity, we assume the tunnel couplings between the 
molecular QD and the two leads to be symmetric F^ = 
Ffl = r. Therefore, one can rewrite the resulting Eq. (19) 
in a compact matrix form: 



[(u + iT)! ~ B]G''iuj) 



(21) 



in which B is a N x N symmetrical tridiagonal matrix 
withi3„„ = ed+nhwph, i?„(„-i) = -X^/n, andi?„(„+i) = 
— Avr^ + 1. A similar result is obtained for the advanced 
GF, Ciuj), and we have G° [C]'' in w-space. The 
element G„„ corresponds to the probability amplitude of 
propagating in a state that starts and ends with the same 
number of n phonons, while an element G,„„ (n ^ m) 
denotes the probability amplitude of starting with m and 
ending with n phonons. 

Furthermore, employing the Fisher-Lee relation:^^ 



s{uj) 



(22) 



we immediately obtain the transmission and reflection 
amplitudes in the presence of a phonon bath: 



SL(R)m,L{R)n — ^ L(R)mn — ^mn 



-Srnn + irG^„„(w), 
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SL{R)m,R{L)n — tR{L)mn ~ Unn ~ *rG^„„(w). 



(23) 



It is easy to check that they obey properties (1) r^r + 
t^t = /; (2) rH + t^r ~ 0, due to the unitary property 
of the scattering matrix; and (3) r = —I + t; (4) 2tH = 
+t. 

Inserting the resuhing transmission and reflection am- 
pHtudes, Eq. (23), into Eqs. (12) and (15), we obtain, 
after some algebra: 

SlUO) - Srr{0) = ^ / dujJ2{[it^t)mn? 

mn 

X [P„/£(a;)(l - /r(^)) + P™/]?(c^)(l - /]?»)] 

+P™/]5'(a;)(l-/r(c.))]}, (24a) 

and Slr[Q) = Srl{Q) = — S'll(O). We can also deter- 
mine the elastic contribution to shot noise by imposing 
the constraint of elastic tunneling m = n: 



^ll(O) 



2e2 



xR 



h I {(*^*)n« 

n 

P4/£(c.)(l - f^W)) + /]J(a;)(l - /£(c.))] 
[ithUrPnm^)-fRM'}- (24b) 



Thus, the total current II and its elastic part are given 

by.21,24 



X {P„/r(c^)[l - f^i^)] - P™/]?(c.)[l - m^^)]} ,(25a) 
and 



y d^^(itt)„„p„[/£(u;)-/]^(c.)]. (25b) 



Furthermore, the Fano factor, which measures the devi- 
ation from imcorrelated Poissonian noise, is defined as: 



F = 



Sll{0) 
2elL 



(26) 



B. Numerical results and discussions 



In this subsection we address the numerical calculation 
of zero-frequency shot noise through a single-molecular 
QD in the presence of a phonon bath, on the basis of 
Eqs. (21), (23) and (24). For simplicity, we assume the 
bias voltage is distributed symmetrically, i.e., +V/2 on 
the left lead and —V/2 on the right lead. We also set the 
phonon energy Uph = 1 as the unit of energy throughout 
the rest of the paper and choose the Fermi levels of the 
two leads as the reference of energy /xl = fJ.R = at 
equilibrium. Other parameters in our calculations for 
the single-molecular QD are: A ~ 0.5, F = 0.04. and 
T = 0.04. 



It should be emphasized that the present approach is 
an exact method to deal with an interacting electron- 
phonon system having arbitrary coupling strength if the 
maximum number of phonons Nph involved in the cal- 
culation is taken up to infinity. Unfortunately this is 
impossible in a real calculation and we have to truncate 
this number to a finite value chosen to insure computa- 
tional convergence with desired accuracy. The appropri- 
ate value of Npfi depends on the energy of the Einstein- 
phonon mode, the EPI constant, and the temperature 
of the system under investigation. For the parameters 
involved in the present paper, we choose Nph = 8 to ob- 
tain results with accuracy better than 1% in this and 
next sections. 

In Fig. 2 we plot the total and elastic currents (a) and 
shot noises (b) as functions of the bias voltage V. For 
comparison, we also plot the current and shot noise for 
the same system without the EPI as dashed lines. Obvi- 
ously, the shot noise and current have similar step char- 
acteristics with increasing bias-voltage. When the elec- 
tron is not coupled to the phonon mode, only one step 
structure occurs when the bias- voltage matches the Fermi 
energy of the left lead with the single level in the QD, 
leading to the single peak in the differential conductance 
dl/dV and the differential shot noise dS/dV shown in 
Fig. 2(c). In the presence of a phonon bath, we find 
that: (1) The overall spectra are both shifted by \^ /ujph 
due to the polaron effect (for instance, the main peaks 
in dl/dV vs. V and dS/dV vs. V corresponding to reso- 
nance with the QD level shift from the position V/2 = ed 
to td~ l^ph, with slightly suppressed amplitudes); and 
(2) new resonant peaks separated by the frequency of the 
phonon mode, hujph, appear due to emission of phonons, 
indicating that a new pseudochanncl has opened and par- 
ticipates in contributing to tunneling. 

Figure 3 shows the bias-dependent Fano factor for the 
system under study. Roughly speaking, it is known 
that the current fiuctuation for a noninteracting two- 
terminal conductor is ^ ^'0 tempera- 
ture, and the corresponding Fano factor F is proportional 
to X]n^n(l ~ ^n)/Sn-^"' whcrc T„ is thc transmission 
probability of channel n.^^'^^ Therefore, there is no cur- 
rent noise if the transmission probability is r„ = or 1, 
and the Fano factor is correspondingly equal to P = 1 
or 0, respectively. This is case in the low bias- voltage re- 
gion, below the single-level resonance. In contrast, in the 
high bias- voltage limit, the Fano factor tends to ^ for a 
symmetric noninteracting system, T l = Ffl. In addition, 
we observe that the Fano factor is significantly enhanced 
and exhibits more fine structure due to phonon effects in 
this region, in comparison with the uncorrelated noise. 



IV. DOUBLE QUANTUM DOTS 

Finally, we examine nonlinear transport and fluctua- 
tions in semiconductor CQDs in a series arrangement. 



for which 



V^^ = 0. As shown in Fig. 1(b), the 
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FIG. 2: (a) The calculated total current I (thick line), elas- 
tic current (thin line); (b) the zero-frequency shot noise 
S (thick line), its elastic part S"^' (thin line); and (c) the cor- 
responding diflterential conductance dl/dV (dashed line) and 
differential shot noise dS/dV (solid line) as functions of ap- 
plied bias voltage for a single-molecular QD with bare level 
ed ~ 1.0. For comparison, we also plot the current and shot 
noise in the case without the phonon bath. The other pa- 
rameters used in calculation are: F — 0.04, ujph — 1.0, and 
A = 0.5. The temperature is set to T = 0.04. 

Hamiltonian for the two-site case is given by: 

Hd^Y. ^od]d,j - Tc{d\d2 + H.C.), (27) 

where denotes the tunnel coupUng between the two 
QDs. In terms of the joint electron-phonon representa- 
tion, it becomes: 

Hd = ^ (ej n)(j,n| 

3,n>0 

-^Te(|l,n)(2,n|+H.c.). (28) 

ri>0 




0.4 I ' 1 ' 1 ' 1 ' 1 ' 1 

1 2 3 4 5 

V 

FIG. 3: Fano factor vs. bias-voltage V . Solid line denotes 
results with EPI, while the dashed line denotes results without 
EPI. The parameters are the same as in Fig. 2. 

We also ignore higher-order tunneling processes by as- 
suming that coherent electron transfer between the two 
QDs, governed by T^, does not explicitly involve electron- 
phonon exchange. Correspondingly, the retarded GFs 
can be written in a block-matrix form: 

/Gil (a;) Gl^{u)\_ 

({^ + ^T,)I-B, C 

\ C iiu + iTR)I-B2 J ' ^^^^ 

in which Bj {j ^ 1,2) is a N x N symmetrical tridiag- 
onal matrix: Bj^nn = fj + nJiiOph, -Bj,n(n-i) = -X-j^/n, 
and „(„_|_]^) = —Xj\/n + 1, and C is a iV-dimensional 
diagonal matrix Cmn — —SmnTc- Once again, we can 
write the transmission and reflection amplitudes through 
the CQDs in terms of the retarded GFs, G'y(tj): 

SL{R)m,L{R)n — rL{R)nm — —^mn + *rG'5^;^(22),mn('^)' 
SL{R)m,R{L)^i ^ ^R{L]mn = *rGi2(21) (w) . (30) 

Considering the properties of the matrices involved, we 
simplify the total current II as: 

mn 

- /EM)} 

= if d0jY.(^W),nn{Pnm^){l~rA'i^)) 

mil 

-P™/]J'(c.)(l-/2(c.))}. (31) 

Furthermore, the zero-frequency shot noise S'll(O) is 
given by 

5ll(0) = ^ / duY,[\{AHUn? [Pnm^) 
run 

x(i - friLo)) + p^rj^{Lo)a ~ /]j(c.))] 

+ |(4rL)mnr[Pn/rH(l-/]?M) 

+P™/]JXc.)(l-/£(c.))]}, (32) 
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and also Slr{0) = 5hl(0) = -^ll(O) = -5flfl(0). 
Their elastic contributions can be obtained by setting 
m — n. 

To compare this theoretical analysis with experimental 
results, we numerically calculate the current and shot 
noise based on Eqs. (31) and (32) as functions of the 
energy difference between the two QDs, e = ei — €2, which 
can be tuned by varying the gate voltage. We address 
only the case of high bias-voltage, V = fiL — fJ-R ^ 0. 
While the present formalism is valid for arbitrary choice 
of coupling parameters, tunneling strengths T and Tj,, as 
well as EPI constants Ai(2) , wc specialize in the interest of 
comparison with the CQD experiments,^^ choosing these 
parameters as: Tl = Tr = T = 0.25, = 0.25, Ai = 
— A2 = A = 0.1875, and temperature T = 0.125 {ujph = 
1). In our calculations we set eV = 10 so that hl 
ei,e2 3> /ifl, in which case nonlinear transport does not 
significantly depend on /zl and jiR at low temperatures. 

In Fig. 4(a), we plot the resulting currents, including 
the total current, its elastic and inelastic components, 
and the current without EPI, as functions of energy dif- 
ference, £. The external voltage, V , gives rise to a finite 
tunneling current through the CQDs. In the case without 
EPI, one can observe a single maximum in the current 
spectrum at e = 0, i.e., the resonance point at which 
the energy level of the left dot is matched with that of 
the right dot, and there is a rapid decrease upon depar- 
ture from this point. In contrast, the inelastic current 
due to phonon coupling is remarkably enhanced on the 
positive branch of the energy difference, e > 0, in which 
phonons are emitted during the tunneling process. One 
can even unambiguously observe some peaks when the 
condition for spontaneous phonon-emission-assisted res- 
onance, e = nhuph {n = 1,2 in our figure), is satisfied. 
However, for negative energy gap, e < 0, the contribu- 
tion of phonon absorption processes to inelastic current 
is negligible because there are not enough phonons avail- 
able to be absorbed at low temperature. Moreover, our 
calculations predict a shouldcrlikc structure for e < 0, 
which arises solely from elastic channels in the presence 
of the phonon bath. 

It should be mentioned that the asymmetric and oscil- 
latory shoulder behaviors in the current spectrum have 
also been obtained by means of the master equation^'^ 
and the nonperturbative real-time renormalization-group 
(RTRG) method. "'^^ Nevertheless, the assumption taken 
in the master equation method (i.e., weak tunnel cou- 
pling) limits its validity for analysis of the current spec- 
trum at stronger values of interdot coupling, Tc- Albeit 
the RTRG approach has not such limitation, this prob- 
lem was not actually studied in Ref. [14]. Once again, we 
emphasize that our studies also assume the phonon bath 
to be in thermal equilibrium perpetually. Treatment of 
nonequilibrium phonon effects for CQDs is beyond the 
scope of the present paper. However, one of the advan- 
tages of the present scheme is that it treats the internal 
tunneling, on an equal footing with the couplings 
to the leads as well as EPI. Thus, it allows us to calcu- 
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FIG. 4; Predicted tunneling current as a function of the en- 
ergy difference, e. (a) The total current /tot, its elastic com- 
ponent /cl, its inelastic part /inch and the current without 
phonon bath (thin line) vs. e (with F = Tc = 0.25, EPI 
constants Ai — — A2 = 0.1875, and temperature T = 0.125); 
(b) The total currents (thick Unes) and currents without EPI 
(thin lines) for various couplings, Tc, between the two QDs. 



late the current spectrum for arbitrary internal tunneling 
strength, Tc- Such evaluations are exhibited in Fig. 4(b). 
With increasing Tc, the inelastic-resonance- induced os- 
cillating shoulder structure is gradually smeared out, in 
agreement with experimental observations.^^ 




0.01 



FIG. 5: Zero-frequency shot noise vs. e. (a) Total shot noise, 
Stot; its elastic part, Sei; its inelastic part, Sinci; and shot 
noise with no EPI (thin line). The parameters are the same 
as in Fig. 4(a); (b) Shot noise with (thick lines) and without 
(thin lines) EPI for several different values of Tc. 

In regard to zero-frequency shot noise, our evaluations 
are summarized in Fig. 5. In Fig. 5(a), where we plot the 
calculated shot noise spectra for the same system studied 
in Fig. 4(a), we see that the shot noise without EPI (thin 



9 



curve) exhibits two peaks located symmetrically around 
the resonant point, e = 0. This behavior is due to the 
fact that there is no noise generation when the transmis- 
sion probability is Ttr = or 1 , while maximal generation 
occurs between these values. In contrast to this, we find 
that phonon coupling enhances the shot noise peak lo- 
cated at the positive e-side (e > 0), but reduces the peak 
on the negative side, leading to breakdown of the sym- 
metric structure. Moreover, interference of the EPI in 
the CQDs gives rise to oscillatory shoulders. It should 
be noted that, since the transmission probability corre- 
sponding to these inelastic resonances is small, the shot 
noise spectrum exhibits only a single peak rather than 
two peaks in the phonon-assisted resonant tunneling re- 
gion. 
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FIG. 6: Fano factor vs. e. (a) Fano factors with (thick 
curve) and without (thin curve) EPI for the same system as 
in Fig. 4(a); (b) Fano factors for several different values of Tc. 



The resulting Fano factor is plotted as a function of 
the energy gap, e, in Fig. 6. At the clastic resonance 
point, e = 0, the smallest Fano factor is observed, im- 
plying that quantum coherence significantly suppresses 
shot noise. For inelastic resonant tunneling, spontaneous 
phonon emission processes also reduce shot noise, which 
induces additional dips in the Fano factor for e > 0. In- 
terestingly, our results predict fine structure in the neg- 
ative branch of e. This is due to the shoulders in the 
current and shot noise spectra stemming from the con- 
tributions of clastic channels. It indicates that the Fano 
factor is quite sensitive to the interaction with phonons. 
Moreover, for large |e|, the shot noise is dominated by the 
Poisson process, i.e., Fano factor F ^ 1. Finally, increas- 
ing the interdot tunneling. T^,, progressively smooths out 
the shot noise spectra and the Fano factors as displayed 
in Figs. 5(b) and 6(b). 



V. CONCLUSIONS 

In summary, we have formulated a method that is well 
suited for the analysis and numerical evaluation of in- 
elastic tunneling and the determination of the nonequi- 
librium zero-frequency shot noise spectra of mesoscopic 
devices in the presence of a dispersionless phonon bath at 
low temperature. Following the original work of Bonca 
and Trugman, we transformed the many-body EPI prob- 
lem onto a one-body scattering problem, by projecting 
the original Hamiltonian in the representation of the 
joint coupled electron- phonon Fock space, i.e., the direct- 
product states of electron states and phonon number 
states. Based on this noninteracting Hamiltonian, we ap- 
plied the Biittiker scattering theory of noise correlations 
in a two-terminal mesoscopic conductor to the case of 
pseudo- multi-channel leads with differing weight factors. 
Furthermore, we derived the formula for zero-frequency 
current-current fiuctuations in terms of transmission and 
refiection amplitudes by properly taking account of the 
appropriate weight factors for the various starting chan- 
nels. 

It is appropriate to remark that the present approach 
to EPI is based entirely on an approximate mapping 
to a one-electron picture providing that the bath is not 
strongly perturbed (if that should occur, then bath-dot 
many-body effects would enter the problem and invali- 
date the one-clcctron picture) . It is believed that this ap- 
proximation is valid only in the high bias- voltage region. 
Indeed, the main purpose of this paper is to analyze zero- 
frequency shot noise for inelastic tunneling under nonlin- 
ear transport conditions. An important advantage of this 
technique is that it does not involve any restrictions on 
the parameters, so that it provides a reliable solution for 
arbitrary EPI constants and timneling strengths to the 
leads, as well as arbitrary internal tunnel couplings. 

Moreover, we have employed the derived formulae to 
analyze the shot noise properties of a single-molecular 
QD and for semiconductor CQDs in a series arrangement 
at low temperature. For this purpose, the NGF method 
provided a convenient tool to calculate the transmission 
and reflection amplitudes in the wide band approxima- 
tion with the help of the Fisher-Lee relation. Some of 
the important findings of our work follow: For the single- 
molecule QD, the differential shot noise displays a peak- 
structure similar to that of the differential conductance 
as a function of increasing bias- voltage due to elastic and 
inelastic resonances. The Fano factor is enhanced due to 
EPI effects in the moderately high bias-voltage region. 
In the case of coupled QDs, both the current and shot 
noise spectra exhibit oscillatory shoulders as functions of 
the energy difference between the two QDs in the positive 
branch, stemming from spontaneous emission of phonons 
in nonlinear transport. In the elastic resonant tunneling 
region, the shot noise exhibits a double peak, however 
only a single peak shows up in the inelastic resonant tun- 
neling region. The Fano factor is quite sensitive to the 
EPI. 
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